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Abstract. We use the stochastic limit method to study long time quantum dynamics 
of a test particle interacting with a dilute Bose gas. The case of arbitrary form-factors 
and an arbitrary, not necessarily equilibrium, quasifree low density state of the Bose 
gas is considered. Starting from microscopic dynamics we derive in the low density 
limit a quantum white noise equation for the evolution operator. This equation is 
equivalent to a quantum stochastic equation driven by a quantum Poisson process with 
intensity 5—1, where S is the one-particle S matrix. The novelty of our approach is 
that the equations are derived directly in terms of correlators, without use of a Fock- 
antiFock (or Gel'fand-Naimark-Segal) representation. Advantages of our approach are 
the simplicity of derivation of the limiting equation and that the algebra of the master 
fields and the Ito table do not depend on the initial state of the Bose gas. The notion 
of a causal state is introduced. We construct master fields (white noise and number 
operators) describing the dynamics in the low density limit and prove the convergence 
of chronological (causal) correlators of the field operators to correlators of the master 
fields in the causal state. 
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1. INTRODUCTION 

The fundamental equations in quantum theory are the Heisenberg and Schrodinger 
equations. However, it is a very difficult problem to solve explicitly these equations for 
realistic physical models and one uses various approximations or limiting procedures 
such as weak coupling, low density, and hydrodynamical limits. These scaling limits 
describe the long time behavior of physical systems in different physical regimes. 

One of the powerful methods to study the long time behavior in quantum theory is 
the stochastic limit method developed by Accardi, Lu and Volovich P . Many interesting 
physical models have been investigated by using this method. In particular, it has 
been applied to study the long time quantum dynamics of a system interacting with a 
reservoir in the case of a weak interaction between the system and reservoir, i.e. in the 
weak coupling limit. It was applied to study the spin-boson model |2j, polaron model 
and nonrelativistic quantum electrodynamics [13 E] , quantum Hall effect [3] , relations 
between Hepp-Lieb and Alli-Sewell laser models [B], bifurcation phenomenon in a spin 
relaxation [7j, etc. 

An important problem is to study the long time dynamics of a quantum system 
interacting with a reservoir in the case the interaction is not weak but the density of 
particles of the reservoir is small, i.e. in the low density limit. To describe a quantum 
physical model to which the low density limit can be applied let us consider an iV-level 
atom (test particle) immersed in a free gas whose molecules can collide with the atom; 
the gas is supposed to be very dilute. Then the reduced time evolution for the atom will 
be Markovian, since the characteristic time t$ for appreciable action of the surroundings 
on the atom (time between collisions) is much larger than the characteristic time £_r 
for relaxation of correlations in the surroundings. The dynamics of the A-level atom 
interacting with the free gas should converge, in the low density limit, to the solution 
of a quantum stochastic differential equation driven by quantum Poisson noise. The 
quantum Poisson process, introduced by Hudson and Parthasarathy jH] (for a description 
of the quantum Poisson process see also Kumerrer [0]), should arise naturally in the low 
density limit, as conjectured by Frigerio and Maassen |10J and later by Frigerio and 
Alicki [TT] . For a general survey of quantum stochastic calculus we refer to the review 
by AttaljE]. 

The quantum stochastic equation for the low density limit was derived by Accardi 
and Lu lEj-jlHj using perturbation series for the evolution operator. A nonperturbative 
white noise approach for the investigation of dynamics in the low density limit is 
developed in [13 Ej, where the mathematical procedure, the so called stochastic golden 
rule for the low density limit, was formulated. This derivation uses the white noise 
technique developed for the case of weak coupling limit by Accardi, Lu and Volovich [T]. 
The approach to derivation of the stochastic equations in [EH [E! is based on use of 
the Fock-antiFock (Gel'fand-Naimark- Segal, or GNS) representation for the canonical 
commutation relations (CCR) algebra of the Bose gas. The approach of the present 
paper does not use the Fock-antiFock representation. 
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We study the low density limit for an iV-level atom (test particle) interacting 
with a Bose gas. Starting from microscopic quantum dynamics we derive quantum 
white noise and quantum stochastic differential equations for the limiting evolution 
operator. A useful tool is the energy representation introduced in (THJ El where the 
case of orthogonal formfactors was considered. In the present paper we consider the 
case of arbitrary formfactors and an arbitrary, not necessarily equilibrium, quasifree 
low density state of the reservoir. To each initial low density state of the Bose gas we 
associate in the low density limit a special "state" (which is called a causal state) on 
the limiting master field algebra. We prove the convergence of time-ordered (or causal) 
correlators of the initial Bose field to the correlators of master fields (which are number 
operators constructed from some white noise operators) in these causal states. These 
states are determined by the diagrams which give nontrivial contribution to the limit. 
The leading diagrams can be interpreted as a new statistics arising in the low density 
limit (new statistics arising in the weak coupling limit is discussed in [T| I18j). 

One of the main results of the paper is that the dynamics in the low density limit 
is given by the solution of a quantum white noise equation, which is equivalent to the 
quantum stochastic equation 

dU t = dN t (S - l)U t (1) 

where U t is the evolution operator at time t describing the limiting dynamics, S is 
the one-particle S matrix describing scattering of the test particle on one particle of 
the reservoir, and N t (S — 1) is the quantum Poisson (number, gauge) process with 
intensity S — 1. The equation describes the evolution of the total system+reservoir 
and can be applied, in particular, to the important problem of derivation of the linear 
quantum Boltzmann equation describing the irreversible reduced dynamics of the test 
particle in the low density limit. Such an equation for the reduced density matrix can 
be easily obtained from the quantum Langevin equation, which can be derived by using 
the quantum stochastic differential equation and quantum Ito table (see sect. EJ) for 
stochastic differential dN t (for a derivation of the quantum Langevin equation see |T7|). 
However, in the present paper we are mainly concentrated on further understanding 
in what sense the Poisson process is an approximation of the usual quantum field 
(Theorem and in mechanism through which the quantum stochastic equation arises 
limit of the usual Hamiltonian equation. 

In order to describe the objects appearing in ((TJ) let us introduce two Hilbert spaces 
TCs and Tl, which are called in this context the system and one-particle reservoir Hilbert 
spaces, and the Fock space T(L 2 {R + ;Tl)) over the Hilbert space of square-integrable 
measurable vector-valued functions from R + = [0, oo) to Tl. With these notations the 
solution of the equation is a family of operators U t ; t > in Hs <8> r(L 2 (R + ; Tl)) (adapted 
process); S is a unitary operator in TCs <%> TC, which is explicitly defined in section 

Let us introduce the notion of a Poisson process. Let X be a self-adjoint operator in 
a Hilbert space JC and ^f(f) the normalized coherent vector in the Fock space T(/C) with 
test function / G K. The number operator is the generator of one-parameter unitary 
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group T(e* AX ) characterized by 

T(e lXX )^(f) = ^{e iXX f); XeR 
The number operator is characterized by the property 

{9(f),N(X)9{g)) = (f,Xg)(y(f),*(g)) 

The definition of N(X) is extended by complex linearity to any bounded operator X on 
/C. Let us consider K of the form L 2 {R + ; TC) = L 2 (M + ) ® 7i. For any bounded operators 
X e B(H S ),X 1 e B(H), and for any t > we define N^Xo^X^ := X Q ®N(x[o,t] 
extend this definition by linearity to any bounded operator K in Tis <S> and call the 
family N t (K); t > of operators in TCs <8> r(L 2 (R + ; 7i)) as quantum Poisson process with 
intensity K. The existence and uniqueness of the solution of the equation in this case 
follows from the general theory of quantum stochastic differential equations. Moreover, 
unitarity of S leads to the conclusion that for each t > U t is a unitary operator (see 
Lemma EJ) • 

For the vacuum state of the reservoir (zero density) such an equation was derived 
by Accardi and Lu |14j. In the present paper we derive this equation for an arbitrary 
quasifree initial state of the Bose gas. The main feature of the present paper is that 
the stochastic equations are derived directly in terms of correlators, without use of a 
Fock-antiFock (or GNS) representation. This simplifies the derivation of the limiting 
quantum white noise equation and allows us to express the intensity of the quantum 
Poisson process directly in terms of one-particle ^-matrix. In our approach the limiting 
equation, the algebra of the master fields, and the Ito table do not depend on the initial 
state of the Bose gas. 

We obtain that the dynamics of the compound system in the low density limit is 
described by: 

1) the solution of quantum white noise equation (jHfjjl or equivalently, quantum 
stochastic differential equation in forms (JTJ), (}4"4"j) and 

2) the family of causal states cpi on the algebra of master fields. 

The reduced dynamics of the system (test particle) in the low density limit for 
the model under consideration, with completely different methods, based on a quantum 
Bogoliubov-Born-Green-Kirkwood-Yvon (BBGKY) hierarchy, has been investigated by 
Diimcke 49], where it is proved that, under some conditions, the reduced dynamics is 
given by a quantum Markovian semigroup. 

In the approach of the present paper the reduced dynamics can be easily derived 
from the solution of the limiting quantum stochastic differential equation. Namely, the 
limiting evolution operator U t and the limit state (pi determine the reduced dynamics 
by 

T t (X) = ip L {U+{X®l)U t ), 

where X is any system observable (bounded operator in 7-^g), (pi(-) denotes partial 
expectation, and T t is the limiting semigroup. This equality shows that U t is a 
stochastic dilation of the limiting Markovian semigroup. Using the quantum Ito table 
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for stochastic differential dN t one can derive a quantum Langevin equation for the 
quantity U^{X ®\)U t . Then taking partial expectation one gets an equation for T t (X); 
in particular, one can get the generator of the semigroup. This is a general feature of 
the white noise approach: one at first obtains the equation for the evolution operator 
of the total system and then gets the reduced dynamics of the test particle. Let us note 
that although quantum stochastic equations, which are derived in [T3J HI] , are different 
from they give the same reduced dynamics. 

The low density limit can be applied to the model of a test particle moving through 
an environment of randomly placed, infinitely heavy scatterers (Lorentz gas) (see the 
review of Spohn [20]). In the Boltzmann-Grad limit successive collisions become 
independent and the averaged over the positions of the scatterers the position and 
velocity distribution of the particle converges to the solution of the linear Boltzmann 
equation. An advantage of the stochastic limit method is that it allows us to derive 
equations not only for averaged over reservoir degrees of freedom dynamics of the test 
particle but for the total system+reservoir. For a rigorous treatment of a classical 
Lorentz gas we refer to PHIZES]- The convergence results and derivation of the linear 
Boltzmann equation for a quantum Lorentz gas in the low density and weak coupling 
limits are presented in [2Z1I2H]- The Coulomb gas at low density is considered in |29j . 

The hydrodynamical limit is described by the Euler equation. In [3U| the Euler 
equation for fermions in the hydrodynamical limit is derived under some assumptions. 

Let us describe the plan of the paper. In section 3 we construct the master fields, 
which are number operators acting in some Hilbert space, and the limit causal states on 
the master field's algebra. We prove that the time-ordered (or causal) correlators of the 
free evolution of the initial field converge in the low density limit to the correlators of 
the master field in these causal states. In section 4 the stochastic Schrodinger equation 
which describes the dynamics in the low density limit is derived. In section 5 we 
bring this equation to the causally normally ordered form. This form is convenient 
for study the reduced dynamics of the system. In section 6 the expressions for the one- 
particle S'-matrix and T-operator are given. In section 7 quantum stochastic differential 
equation for the limiting evolution operator is derived. 

2. AN ATOM INTERACTING WITH A DILUTE BOSE GAS 

Let us explain our notations. We consider a quantum system (test particle) interacting 
with a boson reservoir (heat bath). Let Tis be the Hilbert space of the system. For 
example, for an iV-level atom 7ig = C^. The system Hamiltonian H$ is a self-adjoint 
operator in Tis- The reservoir is described by the boson Fock space T(7i) over the 
one particle Hilbert space Tt = L 2 (M d ) (with scalar product (•,•)), where d = 3 in the 
physical case. Moreover, the Hamiltonian of the reservoir is given by := dT(Hi) 
(the second quantization of the one-particle Hamiltonian Hi) and the total Hamiltonian 
H tot of the compound system is given by a self-adjoint operator on the total Hilbert 
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space Hs <8> r(H): 

H tot := H hcc + H int = H s ®l + l®H R + H mt . 

Here H int is the interaction Hamiltonian between the system and reservoir. The one- 
particle Hamiltonian H\ is the operator of multiplication by some real-valued function 
u(k). The interaction Hamiltonian will be assumed to have the following form: 

H mt := i(D <g> A + (g )A( gi ) - D + ® A + (^)A(^o)) 

where D is a bounded operator in W5, D G B(7Ys); A(g n ) and y4 + (g n ), n = 0, 1, 
are annihilation and creation operators, and go,gi E W are formfactors describing the 
interaction of the system with the reservoir. This Hamiltonian describes scattering of 
particles of the Bose gas on the test particle and can be obtained by quantization of the 
classical interaction potential between particles of two different types with an infinite 
number of particles of one type (particles of the gas) and finite number of particles 
of the second type (test particles). This Hamiltonian preserves the particle number of 
the reservoir, and therefore the particles of the reservoir are only scattered on the test 
particle and not created or destroyed. Such a Hamiltonian was considered by Davies [SI] 
in the analysis of the weak coupling limit. 

The initial state of the compound system is supposed to be factorized: 

p = ps ® y?L,e 

Here ps is an arbitrary density matrix of the system and the initial state of the reservoir 
(fL£ is the gauge invariant mean zero Gaussian state, characterized by 

<plM + U)Mji)) = t (a, T^iL f ) (2) 

for each f,g EH. Here £ > is a small positive number and L is a bounded positive 
operator in TC commuting with St (an operator of multiplication by some function L(k)). 
In the case L = e~@ Hl , so that L(k) = e _/3w( -' c \ where (3 > is a positive number, the 
state ipL£ is just the Gibbs state, at inverse temperature j3 and fugacity £, of the free 
evolution. The fugacity £ = e^; p is the chemical potential. 

The dynamics of the total system is determined by the evolution operator which in 
interaction representation has the form: 

U{t) := e ltH hee e -itHtot 

It satisfies the differential equation 
at 

where the quantity Hi nt (t) will be called the evolved interaction and defined as 

flint (*) =e itH *-H isb e- m *«. 
The iterated series for the evolution operator is 

~, * t n -l 

00 „ „ 

U(t) = 1 + ^(-i) n / dh . . . / d^H^t,) . . . H int (t n ) (3) 

n=l n n 



7 



With the notations 

S t := e ltHl , D(t) := e ltHs De~ itHs 
the evolved interaction can be written in the form 

H- mt (t) := i(D{t) <g) A+(5^o)^(^i) - ® A+^OA^o))- (4) 

We assume the rotating wave approximation 

e itH sDe-' ltHs = D, 

although generalization to the case of arbitrary D is not difficult. 

We study the dynamics generated by the Hamiltonian (J3| in the low density limit: 
n^0,t~l/n(nis the density of particles of the reservoir). The density of particles 
with momentum k in the state tp^^ is equal to 

l-£L(k) 

and goes to zero as £ — > 0. Therefore the limit n — > 0, t ~ 1/n is equivalent to the limit 

£ _> o, t ~ i/e 

Let us consider the time rescaling £ — ► t/£ so that C/(t) — > U(t/£). With the 
notation 

JW*) = ^ + (%/)A(^) (5) 
for any f,g EH., the equation for the evolution operator U(t/£) becomes 

= (D ® A,, f/: i i/) - D + ® N^iWit/Z) (6) 

The reduced dynamics of any test particle's observable X in the low density limit 
is defined by the limit 

hm^(U + (t/0(X®l)U(t/0) 

where ipi,,e{') denotes partial expectation. In [THj it was proved that, under some 
conditions, the limit exists in a small time interval and is equal to T t (X), where 
{It] t > 0} is a quantum Markovian semigroup. The dynamics of the reduced density 
matrix ps(t) is determined through the duality Tr(ps^i(^)) = T^iPsifyX). As was 
mentioned in the Introduction, in the approach of the present paper the limiting 
semigroup can be obtained by using the solution Ut of the quantum stochastic equation 
as 

T t (X) = <p L {U+(X®l)U t ) 

and the generator of the semigroup can be easily derived from quantum Langevin 
equation. The limiting semigroup can be obtained also from quantum Langevin equation 
in [TJj, which is based on a quantum stochastic equation similar to ([!} but much more 
complicated. 

The first step to study the low density limit of the model is to find the limit of the 
field Nf }9 ^(t). This limit we call master fields or number operators. 
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3. THE MASTER FIELDS AND THE LIMIT STATES 

In this section we construct the algebra of the master fields arising in the low density 
limit and the limit causal states on this algebra. We prove (Theorem 1) that time- 
ordered correlators of initial fields (jSJ) converge in the low density limit to correlators 
of number operators constructed from some white noise operators. Theorem 2 states a 
useful factorization property of the limiting causal states. 
It is convenient to use the "projections" 

oo 

P E := -L / dtS t e~ itE = - E) 
lix J 

— oo 

with the properties 

P E P E , = 6(E - E')P El P* E = P E , S t = J dEP E e itE 

For the 5-function of a self-adjoint operator cf. Definition (1.2.1) in pp. 

Let us construct the master space (which is Fock space over some Hilbert space) 
and master fields. For a given Hilbert space Ji and a self-adjoint operator Hi in "H we 
define the Hilbert space X-h,h y as the completion of the quotient of the set 

!^F:R^n s.t. ||F|| 2 :=27r J dE(F(E), P E F{E)) < ooj 

with respect to the zero-norm elements. The inner product in is defined as 

(F, G) = 2n J dE(F(E), P E G{E)). 

We denote by Bf(E,t), B g (E',t') time-energy white noise creation and annihilation 
operators acting in the symmetric Fock space r(L 2 (R + , Xu,Hi)) where L 2 (R + , X^m) 
is the Hilbert space of square integrable functions / : M + — > X nHl . These operators 
(operator-valued distributions) satisfy the canonical commutation relations 

[B g (E,t), B+(E',t>)} = 5(t' - t)5(E> - E)%j(E) (7) 
and causal commutation relations 

[B g {E,t), B+{E',t')] = 5 + (t> - t)5{E> - E) lgJ {E) (8) 
where 5 + (t' — t) is the causal 5-function and 

(g,P E 'f) 



lgJ (E) = J dE', 



i(E> -E-iO) 
%j(E) = 2n(g,P E f) 

In the Appendix we review the definition of the causal 5-function; for a detailed 
discussion of distributions over the simplex and the meaning of two different 
commutators (j2J) and (jHJ) for the same operators we refer to Sect. 7 in |T]. These operators 
are called time-energy quantum white noise due to the presence of 8(t' — t)S(E' — E) 
in (JZj). 
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For any positive bounded operator L in T~L we define the causal gauge-invariant 
mean-zero Gaussian state tpi by the properties (j9*jl- (fT2*Jl : 

for n = 2k VL {B? . ..B?) = J^MB^B^). ..^(Sj B*») (9) 

where the sum is taken over all permutations of the set (1, . . . ,2k) such that i a < j a , 
a = 1, . . . , k, ii < %2 < . . . < ik] := Bf™(E m , t m ) for m = 1, . . . , n, are time-energy 
quantum white noise operators with causal commutation relations (jHJ), and e m means 
either creation or annihilation operator; 

for n = 2k + 1 VLiBl 1 ... B e n n ) = (10) 

MB f (E,t)B g (E',t>)) = <f L (B+(E,t)B+(E',t')) = (11) 

<p L (B+(E,t)B g (E',t')) = x m (t')(9,P E Lf) (12) 

Notice that the "state" (pi does not satisfy the positivity condition. This is a well- 
known situation for the weak coupling limit (see pQ) and is due to the fact that we work 
with time-ordered, or causal correlators. Therefore it is natural to call such "states" 
causal states. 

Definition 1 Causal time-energy white noise is a pair (Bf(E,t),<pi,) , where Bf(E,t) 
satisfy the causal commutation relations (0) and (p^ is a causal gauge-invariant mean- 
zero Gaussian state characterized by Uty- hlty . 

Using the operators B~f(E,t), B g (E,t) we define the number operators as 

N fa (t) = J dEB+(E,t)B g (E,t) (13) 

Finally, for a given Hilbert space 7i and a self-adjoint operator Hi we have the 
following objects: for any £ > the family of operators Nf }9 ^(t) defined by © together 
with the gauge-invariant quasifree mean-zero Gaussian state <^l,£ and the number 
operators Nf )9 (t) together with the causal state (pi. 

The following theorem describes the relation between these objects and states the 
master field in the low density limit. 

Theorem 1 There exists causal time-energy white noise (Bf(E,t),(p£) such thatVn G 
N 

lim^(iV /ljS1 ,e(ii) • ■■N fn>gn ,t(t n )) = <p L {N fugi (h) ■ ■ ■ N fn , 9n (t n )) 

where the equality is understood in the sense of distributions over simplex t\ > t% > 
• • • > t n > 0. The limit causal state ipL is characterized by <f^)-< f7^) and the number 
operators are defined by < TT^) . 

Remark 1 This convergence is called convergence in the sense of time-ordered 
correlators. The fact that we use the distributions over simplex is motivated by iterated 
series (0) for the evolution operator. 
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Proof. Notice that 



where 



N f , g ,dt) = j dEN fa4 (E,t) 
N fM (E,t) := —A+(P E f)A(S m g) 

Therefore 

VL ,d N h, 9 ud t i)--- N U, 9 nAtn)) = [ dE 1 ...dE n <p L4 (N fl ^{E 1 ,t 1 )...N fn)gn£ {E n ,t n )) 



Let us denote for shortness of notation for I — 1, ... ,n, 

At := —j^A + {P E J l ) ] A l := -j=A{S tl/m ) 
In this notation 

r-/.J-V/,,, J /'-I- h)... N fn>gn4 (E n , t n )) = ipL^AfAi . . . A+A n ) (14) 

The state is a gauge-invariant mean-zero Gaussian state. Therefore (|T4*jl is 
equal to the sum of terms of the form 

<PL,t(K A h) ■ ■ ■ ^Mi k A jk )^(A Jk+1 At +i ) . . . ip L> s(A jn Af) (15) 

where k = 1, . . . , n, 1 = i x < i 2 < ■ ■ ■ < ik, jk+i < ■ ■ ■ < jn, H < ji for I = 1, . . . , k and 
ji < ii for I = k + 1, . . . , n. We say that ()15j) corresponds to a nonconnected diagram 
if there exists m G {1, . . . , n} such that i\ < m ji < m. Otherwise we say that (|TKj) 
corresponds to a connected diagram. 

Let us prove that all the connected diagrams except only one corresponding to the 
case k = 1 are equal to zero in the limit. One can write ()15|) as 

1 - t h )Ek + ... + (t in - t jn )E in ]/t) (i k F{E) + 0{e +l ) 

= 1 exp{i[t n (E n - E an ) + ...+ tl (E 1 - E ai )]/C] (^ k F(E) + 0(£» 
= -L exp{z[(t„ - t n ^)u n (E) + ... + (t 2 - tJuiWyt} (e-'FiE) + 0(£ fc )) 

= ^ • • • - — | — (r-^(£) + o(a) (i6) 

where (ai, . . . , a n ) is the permutation of the set (1, . . . , n), uii(E) = E n + . . . + E% — 
E an — ... — E ai for I — 2, . . . , n and 

k n 

F(E) = l[(g n ,P El Lf H ) \{(gmPE H h) 

1=1 l=k+l 

Notice that for a connected diagram all the functions oji{E) are not identically zero. In 
fact, suppose that u m (E) = for some m G {2, . . . , n}. In this case one has the identity 

E m + ... + E n = E am + . . . + E an 



-k+V 
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(where E a , E a > for a ^ a' are independent variables) which means that (a m , . . . , ot n ) is 
a permutation of the set {m, . . . , n} and hence (oti, . . . , a m _i) is a permutation of the 
set {1, . . . , m — 1}. Let us choose any I G {1, . . . , n} and consider the term tj l (Ej l — E it ) 
in the exponent in the second line of JTEJ). If ji < m, then since ii = Oj, and ct^ 
belongs to the set {1, . . . , m — 1} one has %\ = G {1, . . . , m — 1}, and vice versa if 
otj t = %i G {1, . . . , m — 1}, then ji — m — 1. This means that if u)\ are not identically 
zero, then ()15j) corresponds to a connected diagram. 

Let us consider the case k > 1. Then, if (fTK|) corresponds to a connected diagram, 
the functions uji(E) are not identically zero. In this case, since there exists the limit 

lim = 5 + (ti-t l - 1 ) 



and the limit of the product of such terms in (fTo^l . and k — 1 > 0, the limit of (fTBj) is 
equal to zero. 

Now let us consider the case k = 1. In this case (fi3j) has the form 

(f L ^(AtA n )(f L ^(A 1 A+) . ..<f L ^(A n -xA+) = 

- exp{i[(ti - t n )#i + (t 2 - h)E 2 + ... + (t n - t n ^)E n )/^ (tF(E) + 0(£ 2 )) 

{E)H e i(t2-ti)u 2 (E)/Z 



(F(E) + O(0) (17) 



where ooi(E) = E\ — E\. Using the limit (J49)) one finds that the limit of the right-hand 
side (RHS) of (fTTj) is equal to 

5 + (t 2 - tl ) . . . 5 + (t n - t n ^)(g n , P El Lf x ) ^9i,PeJ2) (g„-i,iWn) 
+v ; +v " Ayn ' 1 J 1 i(E 2 - E l - iO) i(E n - E x - iO) 

After integration over . . . E n it becomes equal to 

5 + (t 2 -t 1 )...5+(t ri -V 1 ) y dE{g n ,P E Lf 1 )-y 9ltfa (E)... lgn _ 1 j n (E) (18) 

This proves that only one connected diagram survives in the limit. 
Now let us consider the quantity 

MN h Ah)---N fn , 9n {tn)) 
With the notation 

B+ := B+(Eh,ti); B t :— B gi (E h ti), 
it can be written as 

J dEt . . . dE nVL (B+B 1 . . . B+B n ) (19) 

Notice that on the simplex t± > t 2 > • • • > t n > causal ^-functions 6 + (ti +m — ti) for 
m > 2 are equal to zero. Therefore for m > 2 one has ip L {B tl Bf i+m ) oc 5 + (t l+m — ti)=0 
and hence the integrand in (fTT?|) can be written as 

n-l 

( PL (B+B 1 . . . B+B n ) = M B t B k)M B l B }) ■ ■ ■ M B k-lBi) VL (K+l B k+l ■ ■ ■ B n B n) 

k=l 

+^ L ( J B+ J B n )^ L ( J B lJ B 2 +) . . . <p L {B n _ x B+) (20) 
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The terms in the sum correspond to nonconnected diagrams. The last term corresponds 
to a unique nonzero connected diagram. Moreover 

J dE x ... dE n (p L (B+ B n )(p L (BiB2 ) . . . ip L {B n . 1 B+) = 5 + (t 2 - h) . . . 6+(t n - t n ^) 

x J dE < g n , P E Lf l > j gi j 2 {E) . . 

which is equal to (fTHjl . 

For n — 1 the statement of the theorem is clear. In fact, 

^i Pl ^(N fa4 (t))=^(g, T ^ L f\ = (g,Lf) = J dE^ L {B^{E,t)B g {E,t)) 

Then proof of the theorem follows by induction using the fact that only one connected 
diagram survives in the limit. 

Remark 2 The fact that in each order of iterated series only one connected diagram 
survives in the limit can be interpreted as emergence of a new statistics ( different from 
Bose) in the low density limit. For a discussion of new statistic arising in the weak 
coupling limit we refer to JIJ/ (see also fTftj). 

The following theorem is important for investigation of the limiting white noise 
equation for the evolution operator. 

Theorem 2 The limit state (pi has the following factorization property: Vn £ N, 
<p L {Bf(E,t)N hai (ti) ■ ■ ■ N fn>gn (t n )B g (E,t)) 

= (p L (Bj(E, t)B g (E, t)){p L (N fugi (ti) . . . N fn , gn (t n )) (21) 

where the equality is understood in the sense of distributions over simplex t > t\ > ti > 
... > t n > 0. 

Proof. From Gaussianity of the causal state ifL (property Q) it follows that 
<p L {Bf{E,t)N h M(t x ) . . . N fn>gn (t n )B g (E,t)) 
= <p L (B+(E,t)B g (E,t))p L (N fl , gi (ti) ■ ■ ■ N fn , gn (t n )) 

+ J dE 1 ...dE n Y^MB^E,t)B gt (E i ,U))---MB^(E j ,t j )B g (E,t)) 
The sum is equal to zero since the last multiplier 

P^Bf. (Ej, tj)B g (E, t)) = X[o, tf ](*) < 9,PeM > 
is equal to zero almost everywhere on the simplex t > t\ > t 2 > . . . > t n > and hence 
is equal to zero in the sense of distributions on the simplex. This proves the theorem. 

Theorem 1 allows us to calculate, in particular, the partial expectation of the 
evolution operator and Heisenberg evolution of any system observable in the low density 
limit. In fact, partial expectation of the n-th term of the iterated series for the evolution 
operator (JHJ) (or equivalent series for Heisenberg evolution of a system observable) after 
time rescaling t — > i/£ includes the quantity 

t *n— 1 

dh . . . [ dt^L^Nf^ih) . . . N fntBnA {t n )) 
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(where f a , g a are equal to g or g\). The limit as £ — > of this quantity can be calculated 
using Theorem 1. For example, the contribution of the connected diagram is equal to 

t t\ t2 *n— 1 



n-1. 



J dti J dt 2 5 + {t 2 -t 1 ) J dt 3 5 + (t 3 - t 2 ) . . . J dt n 5+(t n -t 



x J dE(g n ,P E Lf x ) lguh {E)... lgn _ ufn {E) 

= t J dE(g n ,P E Lf 1 ) lgi ,f a (E)... lgn _ ufn (E) 

Similarly one can calculate the contribution of nonconnected diagrams (they give terms 
proportional to higher orders of t). Summing over all orders of the iterated series 
one can find the reduced dynamics of the system. But in the present paper we will 
get the limiting dynamics in a nonperturbative way, without direct summation of the 
iterated series. This procedure includes derivation of the white noise equation for the 
limiting evolution operator and then bringing this equation to the causally normally 
ordered form. After that one can easily find, for example, the reduced dynamics of the 
system. For the weak coupling limit such a procedure was developed in pQ. A nontrivial 
generalization to the low density limit was developed in JSJEI; where the derivation 
is based on the Fock-antiFock representation for the CCR algebra of the Bose field 
determined by the state <Pl,z- The approach of the present paper does not require a 
GNS representation and is different from approach of fOJ HZ] • 

4. THE WHITE NOISE SCHRODINGER EQUATION 

In this section we derive, using the results of previous section, the white noise 
Schrodinger equation for the limiting evolution operator. 

The evolution operator {/(£/£) satisfies equation (jHJ) which can be written as 

du(t/o 



dt 

where 



-&t(t)U(t/t), 



H&) = i(D ® N go , gU (t) -D + ® N gugo ,s(t)) 

The results of the preceding section allow us to write the limit as £ —>■ of the 
Hamiltonian H^(t). In the notation (|T3~j) the limiting Hamiltonian is the following 

operator in H s ® r(L 2 (M + , X HiHl )): 

H(t) = i(D ® N go>91 (t) -D + ® N gi>go (t)) 

= i J dE(D®Bl(E,t)B gi (E,t)-D + ®B+(E,t)B go (E,t)^ (22) 

The dynamics of the total system (system+reservoir) in the low density limit 
£ —>■ is given by a new evolution operator U t which is the solution of the white 
noise Schrodinger equation 

^ = -\H{t)U u U = 1, (23) 
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or equivalent integral equation 

t 

U t = l + J dh(D ® N^fa) - D + ® N^tSjU^. (24) 
o 

5. NORMALLY ORDERED FORM OF THE WHITE NOISE EQUATION 

Our next step is to bring the white noise Schrodinger equation to the causally normally 
ordered form (Theorem 3), i.e., the form in which all annihilation operators are on the 
right side of the evolution operator and all creation operators are on the left side. Such 
a form is convenient for study of the limiting dynamics (see remark 3 and text after 
remark). In particular, it can be used for derivation of (linear) Boltzmann equation. 
We assume that for each E el, the inverse operators 

T (E) := (l + l go , 9l (E)D+ - lgugo (E)D + (j go , 90 j gugi - lgi , golg{hgi ){E)DD + ) ' 

Ti(E) := (l + l gom {E)D+ - lgi , go {E)D + ( lgo , golgugi - lgimlga , gi ){E)D+ d) ' 
exist. 

Lemma 1 If the evolution operator Ut satisfies with H(t) given by h2ty then one 
has 

B go (E,t)U t = 7go , go (E)T (E)DU t B gi (E,t) + T (E)(1 - lgugo (E)D)U t B go (E,t) (25) 
B gi (E,t)U t = -^(EfrWD+UtBtoiE,*) +T 1 (E)(1 + lgo , gi (E)D + )U t B go (E,t) (26) 

Remark 3 Notice that in the RHS of these equalities the annihilation operators Bf(E, t) 
are on the right of the evolution operator. 

Proof. It follows from © and that 

[Bf'{E,t),Nf y g{t\)\ — 5 + (ti — t)jfj(E)B g (E,t) (27) 

Therefore using the integral equation (}2"4*|) for the evolution operator one gets 

B f (E,t)U t = [B f (E,t),U t ] + U t B f (E,t) 
t 

= J dti(D ® [Bf(E, t), N g0;gi (ti)] -D+® [BfiE^N^ih)])^ + U t B f (E,t) 
o 

= (D Vtgo (E)B gi (E,t) - D + Ugi (E)B go (E,t)y t + U t B f (E,t) (28) 
The second equality in (J28|) holds because, due to the time consecutive principle 

[B f (E,t),U tl ]=0 iarti<t. 
In fact, let us consider the quantity 

t t tn-l 

J dhiBfiE^ufr-V] = (-zf- 1 J dh... J dt n [B f (E,t),H(t 2 )...H(t n )} (29) 
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where the n-th term of the iterated series © for U t has the form 



l 'n- 1 

U[ n) : = H)" J dti . . . y dt n #(t x ) . . . H(t n ) 



The commutator [Bf(E, t), H(tk)] proportional to 5 + (tk — t), hence the commutator 
[Bf(E, t), Hit?) . . . H{t n )} is equal to zero on the simplex t > t± > t? . . . > t n > and 
therefore is equal to zero. 

The third equality in ()28|) holds since from (fTTj) and the definition of causal 5- 
function one has 
t 

j dt 1 6+(t 1 -t)B f (E,t 1 )U tl =B f (E,t)U t 
o 

For a detailed discussion of the time consecutive principle and causal ^-function we refer 
tod. 

After the substitution / = g and / = g\ in (f2HJ) one gets 



B go (E,t)U t = (D lgo , go {E)B gi {E,t) - D+ lgo>gi (E)B go (E,t))U t + U t B go (E,t) 
B gi (E,t)U t = (D lgi , go {E)B gi {E,t) -D+ lgu9l (E)B go (E,t))u t + U t B gi (E,t) 
or equivalent ly 

(1 + lgih91 (E)D + )B go (E,t)U t = lgo , go (E)DB gi (E,t)U t + U t B go (E,t) (30) 
(1 + lgugo (E)D)B 9l (E,t)U t = - lgu9l (E)D + B go (E,t)U t + U t B 9l (E,t) (31) 

After left multiplication of both sides of equality (JBUj) by (l + j gi go (E)D) and both sides 
of (USB by j go7go (E)D one gets 

(1 + lgugo (E)D)(l + lg(hgi (E)D+)B go (E, t)U t = l 9o , go {E)D{l + lgugo (E)D)B gi (E, t)U t 
+(l + lgugo (E)D)U t B go (E,t) (32) 

lgo , go (E)D(l + lgugo (E)D)B gi (E,t)U t = - lgo , go (E)DD + lgugi (E)B go (E,t)U t 

+ lgo , go (E)DU t B gi (E,t) (33) 

Now after substitution of expression (jHEJ) into (JH2I) one has 

(l + l ao AE)D + - l 9uao (E)D + (j go , 9o l 9l , gi - lgugolgom )(E)DD + ^B go (E, t)U t 

= lgo>go (E)DU t B gi (E,t) + (l- lgi>go (E)D)U t B go (E,t) (34) 
One can show by similar computations that 

,90 73i,gi 7gi,go79o,9i 

)(E)D + D^B gi (E,t)U t 

= - lgugi (E)D + U t B go (E,t) + (l + lgihgi (E)D + )U t B gi (E,t) (35) 

Now since we suppose that the inverse operators T Q (E) and Ti(E) exist, we can solve 
the above equations (f3~^) and (|33j) with respect to B go (E,t)U t and B gi (E,t)U t . The 
solutions are given by and (ffilj). and that proves the lemma. 
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Denote 

lgugi {E)DT 1 {E)D + 
lgo , go (E)D+T (E)D 
-DT 1 (E)(l + lgo , gi (E)D+) 
D + T (E)(l- lgugo (E)D) 

Theorem 3 The normally ordered form of equation \23\) is 
dU t 



Ri,i( E ) 

Ro,i(E) 



n,m=U,l 



E / dER m , n (E)B+jE,t)U t B 9n (E,t) (36) 



Proof. Using ()22|) white noise Schrodinger equation ()23|) can be rewritten in a more 
detailed form 

IT T p 

= J dE(D®B+(E,t)B gi (E,t)-D + ®B+(E,t)B go (E,t))U t (37) 

It follows from Lemma ^ that 

D + B go (E,t)U t = R 1A (E)U t B gi (E,t)+R li0 (E)U t B go (E,t) 
DB gi (E,t)U t = —Ro t o (E) U t B go (E, t) - R ,i(E)U t B gi (E,t) 

The statement of the theorem is obtained after substitution of these expressions in (IHTj) . 

Remark 4 An immediate consequence of Theorem 2 is the following factorization 
property of the limiting state ip^: 

<p L (B+jE,t)U t B gn (E,t)) = MBt(E,t)B gn (E,t))MUt) 

This property of the state tpi similar to the factorization property of the state determined 
by a coherent vector \\^f\\ = 1: 

(*,B+jE,t)U t B gn (E,t)S>) = (V,B+jE,t)B gn (E,t)S>)(S>,U t *) 

which is usually used to define quantum stochastic differential equations (the general 
notion of adaptedness and adapted domains which are much larger than the coherent 
ones is given in /. 



Taking the partial expectation of both sides of equation (j3Hj) in the state ipi, using 
the factorization property and noticing that 

<p L (B+jE,t)B gn (E,t)) = (g n ,P E Lg m ), 

one gets the equation 

^ = -TMU,), (38) 
where T is being called drift and is equal to 



£ f dER m>n {E)(g n ,P E Lg m ) 

™ — n 1 *J 



n,m=0,l 
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The solution of (|38j) is 

^L(^) = e- r * 

In the case of orthogonal test functions, i.e. {go, Stgi) = this expectation value for 
the evolution operator was obtained in J^l- Let us note that the expectation value is 
obtained in a nonpertrubative way, without direct summation of the iterated series for 
the evolution operator, and is a result of the procedure of causal normal ordering. 



6. ONE-PARTICLE T OPERATOR AND S MATRIX 

In the low density limit the role of multiparticle collisions is negligible and the dynamics 
of the test particle should be determined by the interaction of the test particle with one 
particle of the reservoir. In the present section we give the expressions for the one- 
particle T-operator and S'-matrix. In the next section we will rewrite normally ordered 
white noise equation (jHfij) in a form of the quantum stochastic equation (f4*4*|) and show 
(Theorem 5) that the coefficients of this equation can be expressed in terms of the 
one-particle S-matrix. 

Because of number conservation, the closed subspace of Hs ® r(7"0 generated by 
vectors of the form u <g> A + (f)§ (u G Tis, f G Ti. — L 2 (W d ), $ is the vacuum vector), 
which is naturally isomorphic to Tis <S> Tl, is globally invariant under the time evolution 
operator exp[i(iYs ® 1 + 1 ® Hr, + V)t]. The restriction of the time evolution operator 
to this subspace corresponds to the evolution operator on 7ig ® TC given by 

exp[i(F s ® 1 + 1 ® Ht + Vi)t\ 

where 

V 1 =i{D®\g Q ){g 1 \-h.c.) (39) 
The one-particle M0ller wave operators are defined as 
fi± = s - lim exp[i(if s + H x + V x )t] exp[-i(i7 s ® 1 + 1 ® H{)t] 
The one-particle T-operator is defined as 

T = V\Vl + (40) 
and the one-particle S'-matrix as 

S = (41) 

Theorem 4 For the interaction h8ty the one-particle T-operator and S-matrix have the 
form 

T=-i J dER m , n (E)®\g m ){P E g n \ (42) 

n,m£{0,l} 

S = 1 - 27T J2 I dER mA E ) ® \P E g m ){P E gn\ (43) 
n,me{0,l} 

Proof. For the case {go, S t gi) = equality was proved in pT] . The proof of (jlSJ) 
and ()43|) for the general case can be done in a similar way. 

Expression ()43|) will be used in the next section for derivation of equation (|47|) . 
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7. QUANTUM STOCHASTIC EQUATION FOR THE LIMITING 
EVOLUTION OPERATOR 

Normally ordered white noise equation equivalent, through identification 

B+(E,t)U t B n (E,t)dt = 2rcdN t (\P E g m )(P E g n \)U t 
to the quantum stochastic differential equation 




where N t is the quantum Poisson process in r(L 2 (M + ) (g> 7i) defined by N t (X) : = 
N(x[o,t] ® X), if X is an operator in H. The stochastic differential dN t satisfies the 
usual Ito table 

dN t (X)dN t (Y) = dN t (XY), (45) 

where X, Y are operators in 7i, and the limit state tp E characterized by the property 

M^dN t (\P E f)(P E g\)) = (g,P E Lf)dt 

The coefficients of quantum stochastic equation (|4^|) can be expressed in terms of 
one-particle S-matrix describing scattering of the test particle on one particle of the 
reservoir. To show this we will use Hilbert module notation. For any pair of Hilbert 
spaces X , A\, if N t denotes the Poisson process on the Fock space T(L 2 (M + ) Cg> A\), then 
for bounded operators X e B(X ), X\ G B(Xi), the Hilbert module notation is [TO] : 

N t (X ®X 1 ) := X ® JV t (Xi) 

With this notation equation (j4^j) can be written as 




An immediate conclusion from (|4*H|) and (|4T)|) is the following theorem which is one of 
the main results of the paper. 

Theorem 5 The evolution operator in the low density limit satisfies the quantum 
stochastic equation driven by the quantum Poisson process with intensity S — 1: 

dU t = dN t (S-l)U t (47) 

Equation (|4*Tj) describes the dynamics of the compound system in the low density 
limit. Using this equation and the Ito table for stochastic differentials one can obtain 
a quantum Langevin equation for the Heisenberg evolution of any system observable. 
Then the corresponding master equation or, equivalently, quantum (linear) Boltzmann 
equation for reduced density matrix of the system can be obtained simply by taking the 
partial expectation of this Langevin equation in the causal state ip E - 

Lemma 2 The solution of f^7| ) is unitary. 
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Proof. Let us show that d(U^U t ) = 0. The operator satisfies the equation 
dU+ = U+dN t (S + - 1) 

One has 

d{U+U t ) = dU+U t + U+dUt + dU+dUt 

= U+dN t (S + - l)U t + U+dN t (S - l)U t + U+dN t (S + - l)dN t (S - l)U t 
Using the Ito table P5]l one gets 

dN t (S + - l)dN t (S - 1) = dN t {(S+ -1)(S- 1)) 
This and unitarity of S leads to 

d(U+U t ) = U+dN t (S + - 1 + S - 1 + {S + - 1) {S - l))U t = 

Now it follows from the initial condition U t= o = 1 that, for any t > 0, U^U t = 1. The 
proof of UtU^ = 1 can be done in a similar way. 

8. CONCLUSIONS 

In the present paper we consider the dynamics of a test particle (iV-level atom) 
interacting with a dilute Bose gas. It is proved that the dynamics of the total system 
converges in the low density limit to the solution of the quantum stochastic equation 
driven by a quantum Poisson process with intensity S — 1, where 5* is the one-particle 
scattering matrix. The limiting equation is derived in a nonpertrubative way, without 
use of iterated series for the evolution operator. The derivation is based on the white 
noise approach and on the procedure of causal normal ordering developed for the weak 
coupling limit by Accardi, Lu and Volovich P^. The novelty of the present derivation is 
that it does not use the Fock-antiFock (or GNS) representation for the CCR algebra of 
the Bose gas, determined by the state y?L,§- This simplifies the derivation and allows us 
to express the intensity of the Poisson process directly in terms of the one-particle S- 
matrix. The notion of causal states is introduced and the convergence of the correlators 
of the free evolution of the initial number operators to correlators of quantum white noise 
operators in causal states is proved. The causal states satisfy the factorization property 
similar to that satisfied by states determined by coherent vectors. This property is 
crucial for study of the reduced dynamics of the system. 
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9. APPENDIX: CAUSAL ^-FUNCTION 

Let us recall the construction for distributions on the standard simplex (cf. pQ). Define 
C := {<p : R + -> C | <f) = a.e.}, 

G\ := {0 : M + — > C | is bounded and left — continuous at any t > 0}, 

C := linear span of {C UCi}. 

For any a > define — a) as the unique linear extension of the map: 
8+(- - a) : <j) G Cx — > 0(a) 

<$+(• - a) : (j) G C -> 0. 

In jTj the following results are proved. 

Lemma 3 . In the sense of distributions one has the limit 



e 



i(t'-t)E/\ 2 



lim = 2n5(t' - t)6(E) (48) 

A^O A 2 

Lemma 4 . In the sense of distributions over the simplex t > t' > one has the limit 

i(t'-t)E/\ 2 i 

I i 25^^ = ^ t '- t »i(¥3io)- < 49 > 

The last equality means that for any / G C, g G 5(M), one has the limit 



limj dt'J dE^- 2 f(t')g(E) = f(t) lim J dE? J ^ 



o 



i(E - ie) ' 
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